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Abstract 

A concentrated suspension of lamellar colloidal particles (e. g. clay) is 
modelled by considering a single, uniformly charged, finite platelet confined 
with co- and counterions to a Wigner-Seitz (WS) cell. The system is treated 
within Poisson-Boltzmann theory, with appropriate boundary conditions on 
the surface of the WS cell, supposed to account for the confinement effect of 
neighbouring platelets. Expressions are obtained for the free energy, osmotic 
and disjoining pressures and the capacitance in terms of the local electro- 
static potential and the co- and counterion density profiles. Explicit solutions 
of the linearized Poisson-Boltzmann (LPB) equation are obtained for circular 
and square platelets placed at the centre of a cylindrical or parallelepipedic 
cell. The resulting free energy is found to go through a minimum as a func- 
tion of the aspect ratio of the cell, for any given volume (determined by the 
macroscopic concentration of platelets) , platelet surface charge and salt con- 
centration. The optimum aspect ratio is found to be nearly independent of 
the two latter physical parameters. The osmotic and disjoining pressures are 
found to coincide at the free energy minimum, while the total quadrupole mo- 
ment of the electric double-layer formed by the platelet and the surrounding 
co- and counterions vanishes simultaneously. The osmotic equation-of-state is 
calculated for a variety of physical conditions. The limit of vanishing platelet 
concentration is considered in some detail, and the force acting between two 
coaxial platelets is calculated in that limit as a function of their separation. 
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I. INTRODUCTION 



Charge-stabilized colloidal suspensions have been the object of intense theoretical 
scrutiny, starting with the classic work of Derjaguin and Landau, and of Verwey and Over- 
beek (DLVO) ||, on interacting electric double-layers, which took its roots in the even earlier 
work of Gouy and Chapman [[| on double-layers near infinite, uniformly charged planes. 
Following the latter pioneering contributions, a large body of work has been devoted to the 
planar geometry, which is relevant for electric double-layers near macroscopic electrodes [|J , 
or near membranes of biological interest |J. DLVO theory applies to spherical geometry in 
particular, and yields an effective interaction between the electric double-layers surrounding 
spherical colloidal particles, in the form of a screened Coulomb potential, the validity of 
which has been tested both theoretically |6],|7| and experimentally ||. The cylindrical geom- 
etry, appropriate for (infinitely) long, stiff polyelectrolyte chains, has been investigated along 
similar lines, starting with the work of Katchalsky and collaborators |5,TTi]. The common 
framework of much of this work is Poisson-Boltzmann (PB) theory, which is a mean field 



approximation within the density functional theory of inhomogeneous fluids fllT 



Poisson-Boltzmann theory applies equally well to the description of the inhomogeneous 
distributions of co- and counterions around an isolated charged colloidal particle (or polyion) 
suspended in an ionic solution, and to a concentrated suspension of polyions, provided the 
"cage" of neighbouring polyions is modelled by a Wigner-Seitz (WS) cell of appropriate 
geometry, to which each polyion is confined. Given proper boundary conditions, the Wigner- 
Seitz model reduces the initial many-polyion suspension to the much simpler problem of a 
single polyion confined with its associated co- and counterions to a cell of volume equal to 
the volume per polyion of the dispersion. The shape of the WS cell should reflect the shape 
of the polyion: it will be an infinite slab for a membrane |J, an infinite coaxial cylinder in 
the case of a linear polyelectrolyte 0, or a concentric sphere surrounding a spherical colloid 

Except in the case of spheres, the aforementioned models and calculations deal with the 
case of polyions of infinite extension (e. g. an infinite plane in the case of membranes, or an 
infinite line or cylinder in the case of polyelectrolyte chains). In this paper, we examine the 
case of rigid membranes or platelets of finite size. Restriction will be made to infinitely thin, 
uniformly charged circular or square platelets. These may be considered as a reasonable 
model for dispersions of smectite clay particles, like the natural montmorillonite clays JT3 



or the synthetic Laponite clays [ 14| . While the former are irregularly shaped polygonal 



particles, the latter are, to a good approximation, of a circular shape (disc-like particles). 
The thickness of a real clay particle is of the order of 1 nm, which is much less than lateral 
dimensions (Laponite discs have a diameter of typically 25 nm), so that the picture of 
infinitely thin platelets is acceptable. The main difficulty in a statistical description of such 
platelets lies in the considerable anisotropy of the particles and of their associated electric 
double-layers, compared to the much-studied case of spherical colloids. 

The main objective of this paper is to obtain expressions for the density profiles of co- 
and counterions around a circular or square platelet of finite size, confined to WS cells of 
cylindrical or parallelepipedic geometry. The calculations are carried out within linearized 
Poisson-Boltzmann theory, and the calculated profiles and potential distributions will be 
used to evaluate key equilibrium properties of concentrated lamellar suspensions, including 
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the free energy, the stress tensor, the capacitance and the quadrupole moment of a platelet 
and its associated electric double-layer. 

Preliminary accounts of parts of this work have appeared elsewhere [13-171 . 



II. DENSITY FUNCTIONAL AND POISSON-BOLTZMANN THEORY 

It is instructive to repeat the derivation of the familiar Poisson-Boltzmann equation from 
the free energy functional of the inhomogeneous fluid of co- and counterions contained in 
a WS cell, and subjected to the "external" electric field due to the uniform charge density 
on the finite platelet placed at the centre of the WS cell. In the case of smectite clays, the 
surface charge density a = —Ze/H p (where S p is the area of the platelet and — Ze the total 
charge, in multiples of the proton charge e > 0) is negative, so that the counterions are 
positive, while the co-ions are negative; both are assumed to be monovalent. The local co- 
and counterion densities (or density profiles) are denoted by p~(r) and p + (r), where r is a 
position vector pointing inside the WS cell. If N + and N~ are the total numbers of co- and 
counterions inside the WS cell, the p a (r) satisfy the normalization condition: 

J^p a {r)dh = N a ; « = +,- (1) 

where Q is the volume of the WS cell, which is equal to the average volume per platelet in 
the colloidal suspension. Overall charge neutrality requires that: 

N + - N~ - Z = 0. (2) 

The free energy functional J r [p + ,p~] may be split into the usual ideal, Coulombic, external 
and correlational contributions 111 111 : 



with 



J~ — Fid + ^Coul + J~ext + -^corr (3) 



F id [{p a }} = k B T £ / d 3 r p«(r) [ln(A^(r)) -I. 

01=+,- 



2e Jn Jn 
FeAiP a }] = \ J Q dh p c (r) <p v ( T ). 



, p c (r) Pc(r') 



r — r 



In these equations e is the dielectric constant of the solvent (generally water) regarded as a 
continuous medium ("primitive model"); A a is the de Broglie thermal wavelength of ions of 
species a; p c (r) = p + (r) — p~(r) is the local charge density (in unit of e); and <p v {v) is the 
electrostatic potential at r due to the uniformly charged platelet, the centre of which coin- 
cides with the centre of the WS cell, conveniently chosen as the origin. The correlation part, 
J-" corr , is not known explicitly, but may be expressed within the local density approximation 



ll,6fl. J^corr w iH be neglected throughout, an approximation which is reasonable, as long as 



the local concentrations of co- and counterions are not too high, a condition which may not 
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be fulfilled for the counterions in the immediate vicinity of a highly charged platelet. For a 
complete formulation of the electrostatic problem, the form ([|) of the free energy functional 
must be supplemented by specifying the boundary conditions satisfied by the resulting mean 
electrostatic potential, or its gradient, on the surface E of the WS cell. 

The equilibrium density profiles p a (r) are those which minimize the free energy functional 
(H), subject to the constraints (H) i. e. 

fr°(r) =/i - a = + '~ (4) 

where p a is the Lagrange multiplier associated with the constraint ([I]), i. e. the chemical 
potential of species a. With JF corr = 0, the functional derivatives are easily calculated and 
the optimum density profiles are found to be given by the Boltzmann distribution: 

P ± W = exp{=F/3ey>(r)} (5) 

where (3 = l/(k B T) is the inverse temperature in energy units, and ip(r) is the total electro- 
static potential at r: 

^(r) = ^(r) + - / dV/^- (6) 
e Jn |r — r | 

which satisfies Poisson's equation: 

W(r) = -^ P (r)~ ft (r). (7) 

In eq. (0), q v (r) denotes the (surface) charge density of the platelet. For a given q v (r), eqs. 
(H) and (0) form a closed set, which may be reexpressed as an inhomogeneous nonlinear par- 
tial differential equation for the potential usually referred to as the Poisson-Boltzmann 
(PB) equation: 



Aire 



VV(r) + p+exp{-/3ey?(r)} - p Q exp{+(3ecp(r)} = q p (r). (8) 



4-7T 



Neglect of J-" C orr clearly points to the mean-field nature of PB theory. The prefactors p^ 
are equal to the fugacities exp(/3/i±)/A ± in the case of an open system, where the WS cell 
is in equilibrium with an infinite reservoir which fixes the chemical potentials p±. For a 
suspension of fixed ionic composition (canonical ensemble), the p^ are determined from the 
normalization conditions (]l|). Equation (||) must be solved, subject to appropriate boundary 
conditions on the surface E of the WS cell. If the surface is regarded as a boundary between 
WS cells associated with nearest neighbour platelets, it is natural to impose that the normal 
component of the electric field E = — Vy? vanish at each point on E. In practice, the 
homogeneous version of the PB equation (||) is solved for all positions r outside the platelet 
(r ^ E p ), and the usual discontinuity of the normal component of the field upon crossing 
the uniformly charged surface E p is taken into account. 

The PB equation has been solved numerically for circular platelets of finite thickness (i. 
e. coin-like cylinders), carrying a positive edge charge, in the limit of an infinitely dilute 
suspension (Q — > oo) fl8fl . In this case, the boundary E of the WS cell is pushed out to 
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infinity, the electrostatic potential </?(r) can be chosen to vanish when |r| — > oo, and the 
prefactors p^ reduce to the macroscopic co- and counterion number concentrations n~ and 
n + (n + = n~). An analytic solution of the PB equation (|8]), for vanishing or finite platelet 
concentration n = 1/fi, is available only in the limit of an infinite platelet in a WS slab ; 
this geometry reduces the problem to the classic one-dimensional Gouy-Chapman problem 
For finite platelets, analytic solutions can be obtained only upon linearization of the 
PB equation (^). For this purpose, it is convenient to redefine the prefactors p^ such that 

p ± (r)=pfexp{ T (3e[<p(r)-^]} (9) 

where ip* is a reference potential to be specified. The resulting linearized PB equation reads: 

Att 

VV(r)-^[y,(r)- 7o ] = --g p (r) (10) 
where the squared inverse Debye length k 2 = 1/ , and the constant 7 are given by: 

k\ = 4tt4 (p+ + p~) (11a) 



7o = \Pl ~ P ) — + f = + , _ + <P (Hb) 
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and £ fl = j3e /e is the Bjerrum length (£ fl ~ 0.7 nm in water at room temperature). The 
normalization conditions (P now reduce to: 

/\r± 

n ± = — = Pf f [li/fefo,*-?)] (12) 
where is the mean potential in the WS cell: 



y^T)*. us) 



v a 

A particularly simple choice is thus to linearize the local densities around Tp, i. e. to take 
ip* =Tp ||16| , in which case 

It is worthwhile to note that linearized Poisson-Boltzmann (LPB) theory may be derived 
from the free energy functional (f|) (with T COTT = 0) via the variational principle (Q), provided 
the integrand in the "ideal" contribution jF id , is expanded to second order in powers of the 
local densities p a (r) from their mean n a ||. 

Before turning to the presentation of explicit solutions of LPB theory for specific geome- 
tries, we address the problem of expressing and calculating key macroscopic quantities like 
the free energy, the stress tensor and the osmotic pressure, from the density profiles. 
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III. FREE ENERGY AND PRESSURE TENSOR 



The Helmholz free energy F of the electric double-layer around a colloidal particle inside 
a WS cell is the key thermodynamic quantity which must be evaluated with care [p],|5|,p!9| . 
Within mean-field PB theory, the free energy F may, in principle, be calculated by substi- 
tuting the equilibrium density profiles p a (r), determined via the variational principle (f|), 
into the functional (§) (with JF corr = 0). This expression of the free energy may be cast in 
the standard form 



F = U — TS 

where the internal energy and the entropy are given by [|T 

1 



(15) 



U = - {q p {r) + e[p + {r)-p (r)] } y?(r) d 3 r 



f / [V^(r) 

57T Jn 



d 3 r - — 



•>7T 



j y?(r)V^(r) dS 



(16) 



and 



TS = -k B T ^ I P a (r) [hi(A»))-l]d 3 r 
„~r Jn 



-k B T J2 N«\n(p«Al 

+,- 

{e [p+(r) - p~{ 



ip(r) + k b T [p + (r) + p-(r)] } d 3 v. 



(17) 



In eqs. fll6|) and (0), the p a (r) are the equilibrium profiles, and use was made of eqs. 
([?]) and @ in going from the first to the second line in each of these equations. The 
resulting expression for the dimensionless free energy f3F reads, in terms of the dimensionless 
electrostatic potential $(r) = /3e<^(r): 



(3F= £ A^ln(pX)-^/ s $ ( r )V$(r).dS 



+ 



p (r) - p + (r) $(r) 



P 



Itt£ i 



[V$(r)] 2 d 3 r. 



(18) 



Expression ([18]) , valid within the non-linear PB approximation, involves integrations over 
the WS boundary surface (which do not contribute if the boundary condition of vanishing 
normal electric field is adopted) and over the volume of the WS cell, and is hence not very 
tractable. Equivalent expressions for the free energy can be obtained by considering generic 
charging processes [p]fl p!6| , p!9[ . For a fixed cell geometry (volume and shape), the variation 
in free energy due to infinitesimal variations of the potential ($(r) — > $(r) + 5$(r)) and of 
the Bjerrum length (£ B — > £ B + 5£ B ), is of the generic forr 

5 (BF) = I [$ V (5$) - <5$ V$] • dS + pU 

OTTl a Jt, 



81 



+ 



$<5 



a 



d 2 r 



In (p+Al 8N + + In p Q -Ai 



5N_ 



(19) 
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where a = —Ze/Y, p is the surface charge of a platelet per unit area. Expressions of the free 
energy can be derived from eq. (|19|) by considering various real or virtual charging processes. 
In practice, we have calculated F using a constant Debye length charging process, from a 
situation where the platelet is neutral (Z' = 0) to its final charge [Z' = Z). For uncharged 
platelets, iV o + = N~ = N = Qn s + Z/2, where n s is the salt concentration. At any stage 
of the process, N ± = N ± Z'/2, and eq. fll9|) shows that in an infinitesimal step, during 
which the boundary condition of vanishing normal electric field is enforced, the free energy 
changes by: 



a 



(35 (F) - 

Integration along this path yields: 
F(a)-F(a = 0) 

N n kT In 



$6[~)d 2 T+ ln(p+A^ 



5N + + [in (p o -A 3 _) 



6N_ 



(20) 



da' + 



{N f-Z 2 /A 

K) 2 . 



^kT In 



N + Z/2 ' 

K-z/2. 



(21) 



The result (0), valid for a given WS cell geometry, may be generalized to the case of an 
infinitesimal change 5Q of the geometry (shape and/or volume) of the cell. A calculation 
given in the appendix shows that for fixed numbers of co- and counterions (N + and iV~) in 
the cell, the resulting infinitesimal change in free energy reads: 



5F = -k B T £ f p°(r)d 3 r - ± f [Vyf d\ 

J oil On Joil 



(22) 



where it has been assumed that the electric field has no normal component on the WS 
surface. Under these circumstances, the previous expression for SF is compatible with the 
usual definition of the pressure tensor in a charged medium [^0 



n 



P(r) + f {Ef 



I E(r) ® E(r) 

An 



where E(r) 



-V^(r), and for non-interacting ions, 

P(v) = k B T P a (r). 



(23) 



(24) 



The compatibility of eqs. ( 22[) and (^3|) plus (|24]) is easily verified by considering a small 
local volume change 5Q due to the displacement by an amount 5£ of a small area element <5E 
at the surface of the WS cell. The pressure tensor (|23| ) satisfies the mechanical equilibrium 
condition 



v- n (r) = 



(25) 



which, upon substitution of ( p3|) into fl25|), and use of Poisson's equation (0), is equivalent 
to the familiar force balance equation 



VP(r)=e Pc (r)£(r) 



(26) 
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There is no obvious definition of the macroscopic osmotic pressure of the co- and counte- 
rions within the present WS model. An element dS of the surface of the WS cell is subjected 
to the force: 



dF = II • dS = ( k B T Y p a + —E 2 \dS. (27) 
It is hence tempting to define the osmotic pressure as: 

n = fc s T £ + f ^ s (28) 

where the averages are taken over the total surface bounding the WS cell. The same expres- 
sion for II follows from the volume derivative of the free energy, for a particular infinitesimal 
dilation of the WS cell. The latter is chosen such that for each surface element dS of the WS 
cell, the infinitesimal volume element is dQ = <i£ 81, where 81 is a constant displacement 
along the normal to dE. Under these conditions, we deduce from eq. (|2"2]) that: 



5F = -k B T8ey I p a (r)cffi- ^-81 1 (Vy?) 2 dZ, (29) 



so that: 



rfcr^-sN'" 1 (30) 

Similarly, one may express the disjoining pressure Il d (i. e. the pressure to be applied to 
maintain the parallel platelets of a stack at a distance H which coincides with the height 
of the WS cell) by considering the variation of the free energy (^) upon increasing H by 
an infinitesimal amount 5H ; the corresponding increase in the volume of the WS cell is 
8fl = S8H, where 5* is the cross-section of the cell parallel to the platelet (cf. the prismatic 
geometries considered in the following sections). By proceding as in the case of the osmotic 
pressure, one arrives at the required expression 

d(F/S) s £ r 21 S Jt-S , . 

dH~ = ' 3 P 8tt ^ ^ J = ' ^ ^ 

where the average is now taken over the cross-section S, i. e. on the surface of the WS cell 
parallel to the platelet, and the normal to the platelet is chosen along the z axis. 11^ is thus 
related to the mean uniaxal stress along the normal to the platelets. 

Another instructive quantity is the capacitance C of the electric double-layer associated 
with a platelet. For finite platelets, the definition of C is ambiguous. We define the ca- 
pacitance in terms of the difference (Ay?) 5 ~ s = Tp s ' — Tp s between the potential averaged 
over the cross-section S' of the cell containing the platelet (z — 0), and the corresponding 
average on the surface S midway between two platelets in a stack (z = H/2, cf. Fig. 1 for 
the case of cylindrical and parallelepipedic WS cells). The capacitance is thus defined by 

C (Ay?) 5 '" 5 = a. (32) 



8 



The reciprocal of the capacitance defines a length A c = C~ x which characterizes the thickness 
of the electric double-layer. 

The total charge inside a WS cell is zero, and due to space reflection symmetry, the 
electric dipole moment associated with the charge distribution inside the cell vanishes. The 
first a priori non vanishing multipole moment of the charge distribution is the quadrupole 
moment Q = Q zz = —2Q XX = -2Q yy : 



Q 



tot. 



1 

2 J si 



K(r) + 



P w 



x — y ) dr. 



(33) 



In the two following sections, the quantities defined in this section will be calculated within 
cylindrical and parallelepipedic geometries, for disc-shaped and square platelets. 



IV. CYLINDRICAL GEOMETRY 



We consider first the case of disc-shaped platelets, of radius r . In ref. |l5j explicit 
solutions of LPB theory were obtained in terms of infinite series of Legendre or Bessel 
functions, for spherical and cylindrical WS cells. We reexamine the latter geometry in some 
detail here. The cylindrical WS cell is of radius R and of height H = 2h (cf. Fig. 1), so 
that Q = 2n R 2 h. The electrostatic potential within the WS cell is a function (p(r, z) of the 
cylindrical coordinates r and z, which satisfies the boundary conditions 



d(p(r, z) 



Or 
d(p(r, z) 



0. 



dz 



r=R 



=±h 



0. 



(34a) 
(34b) 



Since these conditions, as well as the discontinuity of the electric field E z {z = ± ) = 
— (d(p/dz) z=0 ± across the disc (r < r ), involve only the derivatives of the potential, one 
may assume Tp = without loss of generality (cf. eq. (|13|)). Under these conditions, the 
solution of the linearized PB equation fllCf) may be expanded in a Bessel-Dini series [ D| : 



<f(r,z) 



(35) 



where y n is the n th root of J\{y) = —dJ (y)/dy = 0, and J and J\ are the Bessel functions 
of th and 1 st order {yi = 0). The resulting differential equations for the coefficients A n (z) 
are easily solved, leading to the final result: 



$(r, z) = fie tp(r, z) = fle^o + 



K D b 



cosh[K D (/i — \z\)] 



RJ smh(K D h) 
A n Ji(k n r ) 



b R^ 2 Vn sinh (h/A n ) J$(y v 



cosh 



h — \z\ 



Jo{k n r) 



(36) 



where b = e/(2n£ B a) is the Gouy length, A n = R/ ' \Jy\ + k 2 d R 2 and k n = y n /R. 

For any given macroscopic density of platelets, n, and hence a given volume Q — 1/n of 
the WS cell, only the product R 2 h is fixed. Since R cannot be less than the disc radius r , 
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the aspect ratio h/r < (27mr^) . At the upper limit of this range, i. e. for R = r , and 
with the chosen boundary conditions, the electrostatic problem reduces to that of an infinite 
uniformly charged plane in a WS slab of width H = 2h. With R = r the potential in eq. 
(|36|) indeed reduces to the first two terms, which are independent of the radial coordinate 
r, i. e.\ 



lim $(r, z) 

R^r 



$(z) = (3e% + 



1 cosh[n D (h 



K D b 



sinh(«; D /t) 



(37) 



The above expression is precisely the solution of the one- dimensional LPB equation for an 
infinite uniformly charged plane in a slab, and the familiar result <&(z) = exp(— n D z) / (n D b) 
is recovered when h — > oo. 

The density profiles calculated from the potential (^) via eq. (jffy are sensitive to the 
aspect ratio h/r for a given cell volume Q. The "optimum" ratio is determined by mini- 
mizing the free energy F with respect to this ratio. F is calculated via the constant Debye 
length charging process, resulting in eq. (|2l|). The part of that expression which depends 
on the ratio h/r for a given cell volume is: 



.4 



(p a '(r,z= 0) dh 



da' = 2-7T 



o 







ip° (r) r dr 



da' 



rrrr / tp a (r) r dr 
lo 



(38) 



where (p (r) = f(r,z = 0) is the electrostatic potential on the disc and the last line holds 
within the LPB approximation only, where the potential can always be expressed in the 
form y? CT (r) = a f(n D ,r). Substitution of eq. (|36|) (with z = 0) into eq. ( p8|) yields 



L A = l Qe 1 i( r °Y 1 

Z 2 P67 ° 2K D b\\RJ tanh(K D /i) 



' ^ y 2 n 4{y n ) tanh(VAn) 



(39) 



n=2 



An example of the variation of F with /i/r is shown in Fig. 2a. In this example, as well as 
under all physical conditions that were investigated, F goes through a minimum for a ratio 
h/r , such that the physical requirement that R > r is satisfied. Moreover, the location 
of the minimum turns out to be practically independent of the charge density a carried by 
the disc, and of the salt concentration n s ; in other words, the ratio h/r depends "only" on 
the WS cell volume Q, or equivalently on the macroscopic density of platelets n. For any 
density n, the system selects an optimal ratio h/r . The variation of this ratio with n in 
shown in Fig. 3. All figures correspond to Eqgs — 78 and T = 300 K. 

The potential (p(r, z) and the resulting charge density profile p c (r, z) are highly 
anisotropic functions, as expected from the platelet geometry. This requires that a large 
number of terms (typically 50 to 400) be retained in the expansion ( PB[ ) to ensure adequate 
convergence. A typical example of equipotential lines is shown in Fig. 4. 

If n' s is the salt concentration in a reservoir which is in osmotic equilibrium with the 
colloidal suspension, i. e. with the co- and counterions inside the WS cell, then the salt 
concentration n s inside the latter is, within LPB theory, related to n' s by HlBf : 



n,. 



n' 



1 + 



Z 2 



rr 



Zn 
2r7 



< 1, 



(40) 



10 



which is an expression of the familiar Donnan effect. 

The values of the potential and its gradient on the surface £ of the WS cell may be used to 
compute the pressures II and Il rf , according to equations ( |28| ) and (|31|). For a given WS cell 
volume Q, the variations of II and 11^ with the ratio h/r differ: 11^ decreases monotonously 
as h/r increases, while II goes through a minimum for h/r close to the "optimum" value 
minimizing the free energy (cf Fig. 2). At the minimum of the free energy, the two pressures 
are seen to coincide (II = ILi) . Taking into account the conservation of the overall volume 
5 (R 2 h) = 0, eq. (|22|) can indeed be rewritten as 

*f = s - ^ pi - iy , (4i) 

so that the extremum condition SF = implies II = ILi. The variation of II = ILi with 
the platelet concentration n and with the salt concentration is shown in Figures 5a and 5b 
for two values of the platelet charge (Z = 100 and Z = 200). After a shallow minimum, 
the pressure II is found to increase with n (for fixed n' s in the reservoir), whereas it drops 
rapidly with increasing n' s , towards the value of the osmotic pressure in the reservoir. These 
tendencies are reminiscent of the experimental and numerical results of Dubois and coworkers 
concerning lamellar phases of "infinite" charged bilayers . 



The potential distribution may be used to evaluate the capacity C from equation 
The characteristic double-layer thickness is easily calculated to be: 

X c = K f^Vtanhf^. (42) 



R) \2X 



The variation of A c with platelet and salt concentrations is illustrated in Fig. 6a and 6b, 
for two values of the platelet charge. A c reduces to the Debye length X D when r — > R and 
h — > oo, which corresponds to the limit of a single uniformly charged infinite plane. 

Finally, as illustrated in Fig. 2, the quadrupole moment of the charge distribution inside 
the WS cell vanishes at the "optimum" ratio h/r , which minimizes the free energy. This 
coincidence is systematic and may be related to an exact result of Gruber et al. , provided 



the WS model yields an accurate description of a regular, periodic stacking of the platelets 



T7|. A partial explanation follows from the calculation given at the end of the appendix. 
When the Boltzmann's weights @ are linearized, the kinetic part of the free energy variation 
(cf. eq. ( |76a| )) as a function of the aspect ratio is proportional to the total quadrupole (cf. 
eq. ([79|)). Fig. 2b shows this correspondance but also shows that the electrostatic part of 
the free energy variation, defined from eq. ( [76 b| ) , vanishes for the same aspect ratio as the 
kinetic one. 



V. PARALLELEPIPEDIC GEOMETRY 



Instead of a cylindrical WS cell, we now consider a space-filling parallelepipedic cell of 
dimensions L*L*H in the x, y and z directions respectively (cf Fig. 1). The potential 
is naturally expanded in plane waves, compatible with the periodic boundary conditions, 
which are equivalent to the condition of vanishing normal component of the electric field on 
the surface S bounding the prismatic WS cell: 



11 



<^( r ) =E^( k ) exp{ik-r}, 



(43) 



with 



{n x ,n y ,n z ) E ZZ . 



In terms of Fourier components, eq. (|10D becomes 



47T, 



where 



(fc 2 + ^) ( ^ (k) _ 7o4o) = _^ (k) 



q v ( k ) = ^J n q v( r ) exp{ik-r}d 3 r. 



(44) 



(45) 



(46) 



Consider first the case of a circular platelet (or disc) of radius r , as in section IV; q v (k) 
is easily calculated to be 



q v {k) = -^Ji(k /f r ) 



(47) 



with 



ki 



2tt 



The Fourier components of the electrostatic potential are then determined by substituting 
eq. fl4"Tj) into eq. ((4^), and inverse Fourier transformation leads to the desired result: 

$(r) = (3e ^(r) = (3e % + ^ £ cos ^ x + k vV) 



(n x ,n y ) e . 
1 



k, 



cosh 




2 D + ^) 1/2 (^- 




sinh 







(48) 



The resulting free energy, as calculated from the constant Debye length charging process 
(0) and O, is: 



z 2 ^ e7 ° i?m 



E 



(ji x ,n y ,n z ) eZ 3 



■/i (k//r ) 
k* (k 2 + ^ 



(49) 



Similarly, the quadrupolar moment may be calculated from eq. fl4~8|) to be 







tot. 

22 



16 







disc 



E 



U2 _|_ „2 p i k 2 

1 ft 2 + K D n;c =0 K X + K D 



k.y 



(50) 



The key finding is that, for given platelet and salt concentrations, the changes in free 
energy and quadrupole moment, induced by the change of topology of the WS cell, are 
practically negligible, in spite of the completely different analytical expressions (cf. eq. (|3~8"D 
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and (f49D), as illustrated in Fig. 7. A similar conclusion holds for the osmotic pressure 
calculated from eq. (|28|) or the capacity evaluated from eq. (|32|) . 

We finally consider the case of a square platelet of side l , within the above prismatic 
WS cell. In this case, the surface charge per unit area is a = —Ze/1%, while the Fourier 
transform of the platelet charge density is now 

The potential is once more expanded in plane waves with the result 

$(r) =0e<p(r) = f3e lo + -^ £ sin (Mo/2) sin (k y l /2) ' 



^ v-w, v -y-u, -/ , , 

(%,%)6l J 



cosh 




l + kiy/\h- 




sinh 







cos (^ x + ^ y) p - /2 — T _ — — — — 1 . ( :,2 ) 



As in the case of a disc in a cylindrical WS cell, one can check from eq. ( p2|) that the 
potential goes over to that of a uniformly charged infinite plane when l — > L (expression 

(ED)- 

The corresponding free energy is now: 

/3 1 16 _ sin 2 (W2) sin 2 (fc^ /2) 

o (n I ,n v ,n»)eX 3 h x h y ' ^ D J 

while the quadrupole moment reads 

Q " -'411 ("D-T^TT " ^ £ • (54) 



Qplatelet ^ 1 n z =l K d n x =l 

Explicit calculations based on these formulae lead again to results which are quite close to 
those obtained for circular platelets under similar physical conditions, as shown e. g. in 
Fig. 7. In particular, the "optimum" aspect ratio is characterized by the equality of the 
osmotic and disjoining pressures (an expression similar to eq. pT| ) holds), as illustrated in 
Fig. 8. However, the minima of the free energy and of the osmotic pressure do not coincide 
any more. 



VI. INFINITE DILUTION LIMIT 

The limit of very low platelet concentration n, for a fixed salt concentration n s (or 
equivalently n' s ) may be derived from the expressions obtained in sections IV and V for 
prismatic geometries, by letting the volume Q of the WS cell go to infinity. In this limit, 
7 = 0, and the series fl3"6"|) and fliS]) go over to the following integral representation of the 
reduced potential 

r ,-OC e -\J K2 D + k2 \ z \ 

$(r, z) = dk J (kr) Mkro) (55) 

b Jo . K 2 +k 2 
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Along the z-axis passing through the centre of the disc (r = 0), eq. (|55D reduces to: 



$0 = 0, z) 



e 



(56) 



which goes over to the familiar exponential solution of linearized Gouy-Chapman theory in 
the limit r — > oo (infinite plane). The reduced potential at the centre of the disc takes the 
value 

$( 0) = ^=y = ) = -L{l-exp(- K „ ro )}. (57) 

This expression will be used in the concluding section to evaluate the range of validity of 
LPB theory. 

The quadrupole moment around circular or square platelets vanishes identically in the 
zero concentration limit (n — > 0) as a consequence of the theorem by Gruber et al. . 

Finally consider the problem of determining the force acting between two circular 
platelets (discs) V\ and V2, in the limit of vanishing concentration (n — > 0), and for a 
given salt concentration (in the limit n — > 0, n s — > n' s ). Let ni) and (r 2 , n 2 ) be the cen- 
tre positions and normals of the two discs, and let ip(r) be the total electrostatic potential 
due to the two discs and their electric double-layers; ip(r) vanishes as |r| ^00. In the LPB 
approximation, ip(r) satisfies eq fllCf), with 7,, = 0, and two source terms 

\ 4:71 r 1 

V 2 -<)^(r) = -— [q Vi (v) + q T2 (v)] . (58) 

Since the discontinuity of the electric field upon crossing one of the platelets (characterized 
by the surface charge density q (r)) is independent of the presence of the other platelet, the 
solution ip(r) of eq. Q5BD is just the superposition of the solutions for each platelet separately: 

<p(r) = <p 1 (r)+<p i (T). (59) 

This simple property is a consequence of the linearization, and does not hold within non- 
linear PB theory. The solution of eq. (|58|), and the resulting co- and counterion density 
profiles may now be used to calculate the pressure tensor ( p3|) for given positions and ori- 
entations of the platelets. The force Fj acting on platelet i (=1 or 2) follows then from 

integration of II over the two faces E+ and E~ of the platelet: 



E + .;E". 

Pi" 1 Pi" 1 



n-ds,. (60) 



Since the normal rij = dSj/|dSj| has opposite orientations on the faces S^j and E~ i( the 
kinetic contributions p a (r) k B T to the pressure tensor (p3f), which are continuous across the 
platelet, do not contribute to the force (|5D]). The normal component of the electric field 
E = — V<£>, however, suffers a discontinuity across the uniformly charged platelet, and hence 
contributes to the surface integral in eq. flBTp. 

The total electric field E in the immediate vicinity of the platelet Vi may be decomposed 
into a discontinuous and a continuous part: 
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E = + En = ±— o n, + E 

£ 



2tt 



(<0 



(61) 



where the + and — signs go with the upper (S^) and lower (£ Pii ) faces respectively. Sub- 
sisting eq. (|61|) into (|60[) leads to the desired result: 



(7 



E. (c) d 2 S 



1,2), 



(62) 



where the integration is now over the surface of the platelet. For symmetry reasons, the 
only non-vanishing contribution of E^ to the surface integral in eq. (|62"D is the electric field 
due to the other platelet and its associated electric double-layer. In the case of two coaxial 
parallel discs, the force Fi = — F 2 is along the common axis (chosen to be the z-axis) and 
the result ([55]) may be used to compute the gradient of ip along Oz. The resulting force is 
easily cast in the form: 



F z (d) 



Ana 2 



7T7\, 



T 2 / \ 1 f d 
Ji (x) — exp < 

x r n 



+ it 2 D r*> dx, 



(63) 



where d is the distance between the two discs. For any finite salt concentration (i. e. non- 
vanishing k d ), the decay of F z with d is essentially exponential as illustrated in Fig. 9. In 



the limit of vanishing salt concentration 
we find in that limit 



0), F z decays like a power-law. For d^> r , 



F z (d) 



d>r 



nr 2 a 



2 r 



d 



3rl 
2d 2 



(64) 



Note that the force ( |63|) is repulsive at all distances. 



VII. CONCLUSION 

While most of the existing theoretical literature on suspensions of charged lamellar par- 
ticles or membranes deals with the simpler one-dimensional problem of infinite charged 
planes, we have examined in this paper the case of stacks of finite-size, circular or square 
platelets corresponding, for instance, to swollen clays. The intractable many-platelet prob- 
lem is reduced to the much simpler problem of a single platelet within an electrically neutral 
Wigner-Seitz cell of appropriate volume and topology. The co- and counterion density pro- 
files have been obtained from analytic solutions of linearized Poisson-Boltzmann equation 
with appropriate boundary conditions on the surface of cylindrical and parallelepipedic cells. 
The relevant characteristics of the electric double-layer and the resulting thermodynamic 
properties have been calculated over a wide range of physical conditions. The most notable 
results may be summarized as follows: 

a) For a given cell volume Q, the system selects an optimal size ratio corresponding to 
the minimum free energy. The selected size ratio is practically independent of the platelet 
charge density a, and of the salt concentration, but varies with the macroscopic platelet 
concentration n. 

b) The osmotic presure II and depletion pressure 11^ due to the co- and counterions have 
qualitatively different variations with the aspect ratio for fixed a, clay and salt concentration, 
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but coincide at the optimum aspect ratio which minimizes the free energy. We have a simple 
explanation for this coincidence but not for the observation that the osmotic pressure exhibits 
a minimum at the same aspect ratio as the free energy, at least in the case of circular platelets. 

c) The total quadrupole moment of the charge distribution in the WS cell vanishes at the 
optimal aspect ratio; this observation may be related to an exact property of neutral charge 
distributions in thermodynamic equilibrium HO]. A partial explanation of this correlation 
has been given at the end of section IV. 

The limitations of the present LPB theory must be underlined. Returning to the expres- 
sion ( |57D of the reduced electrostatic potential at the centre of an isolated charged disc, it 
is clear that linearization of the Boltzmann factors in the PB equation (|8[) is justified only 
provided |$(0)| <C 1, i. e. when |6| > X D (which corresponds to the limit of low surface 
charge a or high salt concentration), or r < X D and r < |fe|. This latter condition is rather 
academic for Laponite clay discs (it yields Z < 10). Fig. 10 summarizes the limits of validity 
of the linearized theory for an isolated platelet. Strictly speaking, the above criteria apply 
in the infinite dilution limit n — > 0, and we may expect that they become necessary but 
not sufficient conditions for finite concentrations n [|, so that the shaded area in Fig. 10 
should decrease. Neglecting the clay contribution to the Debye length, the criterion X D < \b\ 
yields Z < 2.3 10~ n ^ 2 r 2 , for a 10 _n molar monovalent salt concentration, with r expressed 
in nanometers. In practice, we carried out calculations for r = 12.5 nm (a typical size of 
Laponite particles); this means that Z must be chosen less than about 400 10~ n//2 . In fact, 
many of our calculations were carried out for Z = 100, in which case, strictly speaking, 
the salt concentration would have to be 10 _1 M or higher for LPB theory to be applicable. 
However, the LPB results should remain, at least qualitatively, valid at lower salt concentra- 
tions. The discussion shows that more realistic calculations, for values of Z typical of clays 
(e. g. Z ~ 10 3 for Laponite), and lower salt concentrations (to ensure stability against floc- 
culation) will require solutions of the full (non-linear) PB equation within the WS geometry. 
Work along these lines is in progress. 

Acknowledgments: Useful conversations with Marjolein Dijkstra, Paul Madden and 
Thierry Biben are gratefully acknowledged. 
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APPENDIX 



Consider the case of a charged platelet V confined in a WS cell with its neutralizing 
counterions and salt. A local change Sfl in volume of the cell (with or without conservation 
of the overall volume) changes the potential and the microion densities, but not the charge 
density of the platelet. The boundary conditions of vanishing normal electric field on the 
surface £ of the cell are enforced during this change. It will be shown that when the total 
numbers of ions N a in the cell are held constant, the free energy F changes according to eq. 
(P^). From eq. flTB|), the change in the internal energy can be written: 

5U=- f (q5p + p5q)d 3 r + - f qpd 3 r, (65) 
2 J n 2 J sn 

where q(r) = q P (r) + e[p + (r) — p~(r)] is the total charge density. It follows from Poisson's 
equation (|7|) and integrations by parts that 

/ q5(pd 3 r= f p5qd 3 r + — I <pV (5<p) dS. (66) 
Jn Jn An Jt, 

For any vector field A(r), the elementary variation of the flux over the surface £ can be 
written: 

5(j>A-ds}=j>5A-dS + jv-A<Pv. (67) 
With A = ip Vy?, the above relation can be used to differentiate the identity 

j> pV (p) ■ dS = 0. (68) 

One finds 

f (pV (8<p) ■ dS = [ qpd 3 r-—[ (Vpf d 3 r, (69) 
Jt, Jsn An Jsn 



e 

An 



so that SU can be cast in the form 
6U 



f p5qd 3 r + ( qpd 3 r-^- [ (V^) 2 d 3 r. (70) 
Jn Jsn on Jsn 

From eq. flTTD and relation (||), the change in entropy reads: 

SS = k B V / p a (r)d 3 r+ [ pqd 3 r + j p5qd 
a ~^_Jsn Jsn Jn 

-k B £ In (pM) \ [ S P a dh + / P a dh 



j3 r 



(71) 



where it was remembered that 



8q p (r) = and / q v (r) <p(r) d 3 r = 0. 

J sn 

Since SN + = SN~ = 0, the expression for the entropy simplifies to: 
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5S = k B y [ p a (r)d 3 r + [ <pqd 3 r+ [ <p8qd 3 r. (72) 

Upon substitution of eqs. (|70|) and ( [72]) into 5F = 5U — T5S, we obtain eq. (p2|) 

^=-fc B T E / P a {v)dh-^f [V V fd 3 r 
„,~^_Jsn Sir Jsn 



n(r)cfr, (73) 
where the following local osmotic pressure was introduced 

U(r) = k B T J2 P°(r) + f (V^) 2 . (74) 

Near the surface E, the relation between II (r) and the pressure tensor (|23|) is 

n(r)=n E -lT-n s , (75) 

so that the osmotic pressure II may be defined as the average of II(r) over the surface E (cf. 
eq. (H)). 

The variation SF can be partitioned into kinetic and electrostatic contributions 

SF eL = - / [Vy?] 2 d 3 r. (76b) 

Within linearized PB theory, we shall finally prove that, when the shape of the cell is 
modified at constant total volume, the resulting variation of the free energy is proportional 
to the total quadrupole moment of the charge distribution inside the cell. 

With the chosen boundary conditions, the total quadrupole defined from eq. ( |33"D is 

Q1T = \ jf <p{r) V (2z 2 - a; 2 - y 2 ) ■ dS. (77) 
Within linearized PB theory, 

<5F kin - = (3e (p+ - p-) I v (t) d 3 r. (78) 

v 7 J oil 

Consider the case of a cylindrical WS cell. Taking into account the conservation of the 
overall volume, it is straightforward to show that 



H- 



,SF 



kin. 



1 



6H 

A similar result holds for a parallelepipedic cell 



5 7o«iQSf (79) 
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FIGURES 



Fig 1. The prismatic cells considered in this paper. Fig la shows a circular platelet 
in a cylindrical cell, whereas Fig lb represents the parallelepipedic cell containing either a 
circular or a square platelet lying in the dashed region S'. 

Fig. 2. Variations of the free energy F, the osmotic pressure II and the disjoining 
pressure n d with the aspect ratio (upper curves, part a)). For illustrative purposes, the free 
energy has been shifted by an arbitrary constant along the vertical axis. Also shown is the 
total normalized quadrupole (lower part, Fig 2b)). The kinetic and electrostatic parts of 
the free energy variation are defined in the appendix (cf eq. ( |76a|) and (|76b[) ). In order to 
check the validity of eq. (79|), we choose a -1 = 7 k 2 Q^ c /2. The normalization pressure 



il is defined with the macroscopic concentrations of co- and counterions inside the cell: 
Tl = k B T(p+ + p;). 

Fig. 3. Variation of the optimal aspect ratio minimizing the free energy with clay 
concentration. The data shown correspond to a circular platelet in a cylindrical cell, with 
n' s = 10" 3 M. 

Fig. 4. Equipotential lines in the plane (k d t, k d z), with arithmetic spacing between the 
isopotentials (/3eA(p = 0.5 between two succesive curves). Here, Z = 100, n = 10~ 5 M and 
n s = 10~ 3 M. For these parameters, the optimal aspect ratio is h/r ~ 2.0. The reduced 
radii of the disc and cylinder are respectively n D r ~ 1.59 and k d R ~ 4.14. The summation 
in eq. (|36"D was truncated after n max = 50. 



Fig. 5a. Clay concentration dependence of the osmotic and disjoining pressures, eval- 
uated at the optimal aspect ratio minimizing the free energy (where II = ILj). The salt 
concentration in the reservoir is held constant (n' s = 10~ 3 M) and defines the normalization 
pressure n reservoir = 2k B Tn' s . 

Fig. 5b. Variation of the osmotic and disjoining pressures with the salt concentration in 
the reservoir (n' s ). For n' s — > 0, Il/Il res ' diverges like l/n' s , while II — > Cst. 



Fig. 6a. The characteristic double-layer thickness (see eq. ([42])) as a function of clay 
concentration, for the optimal aspect ratio minimizing F . The salt concentration in the 
reservoir is n' s = 10 ~ 3 M. 

Fig 6.b. Variation of the reduced double-layer thickness with salt concentration for 
n = 5.10~ 5 M. 

Fig. 7. Free energy and total normalized quadrupole vs. h/r for discs in cylindrical 
and parallelepipedic WS cell (r = 125 A); same quantity as a function of H/l = 2h/l , for 
a square platelet of identical area and surface charge in a parallelepipedic cell (Z = 221 A). 
The vertical lines emphasize the correlation between the minimum of the free energy and 
the vanishing of the total quadrupole. 

Fig. 8. Free energy, osmotic and disjoining pressures as functions of the aspect ratio 
H/l , for a square platelet in a parallelepipedic cell. The data correspond to n — 10~ 4 M, 
n' = 10" 3 M, Z = 100 and l = 250 A. 
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Fig. 9. Variation of the force F z with the distance d between the centres of co- axial discs 
(F = 47T 2 r„<T 2 /e). The four curves correspond to different salt concentrations. 

Fig. 10. Limits of validity for the linearized approximation of PB theory, for a disc of 
radius r in the infinite dilution limit. The dashed line corresponds to the case of an infinite 
plane (in which case r is an arbitrary normalization length). A D denotes the Debye length 
and b is the Gouy length. 
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